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Abstract. The equihbrium state of a turbulent clumpy gas disk is analytically investigated. The disk consists 
of distinct self-gravitating clouds. Gravitational cloud-cloud interactions transfer energy over spatial scales and 
produce a viscosity, which allows mass accretion in the gas disk. Turbulence is assumed to be generated by 
instabilities involving self-gravitation and to be maintained by the energy input from differential rotation and 
mass transfer. Disk parameters, global filling factors, molecular fractions, and star formation rates are derived. 
The application of our model to the Galaxy shows good agreement with observations. They are consistent with the 
scenario where turbulence generated and maintained by gravitation can account for the viscosity in the gas disk 
of spiral galaxies. The role of the galaxy mass for the morphological classification of spiral galaxies is investigated. 
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1. Introduction 

Galactic gas disks are not continuous but clumpy. The 
structure of the interstellar medium (ISM) is usually hi- 
erarchical (Scalo 1985) over length scales of several mag- 
nitudes up to ^100 pc. The clouds are not uniform nor 
isolated and their boundaries are often of fractal nature 
(Elmegreen & Falgarone 1996). Whereas the atomic gas 
(Hi) is mainly in the form of filaments, the molecular 
gas is highly clumped. The largest self-gravitating molec- 
ular clouds (giant molecular clouds = GMC) have sizes of 
~30 pc and masses of ~ IO^Mq (see e.g. Larson 1981). The 
GMCs have a volume filling factor <j)y ~ 10"''^, with the 
ratio of the diameter of a typical GMC to the vertical scale 
height of the GMC distribution (=130 pc, Sanders et al. 
1985) ^0.4. In this respect, the GMC distribution resem- 
bles more a planetary ring (Goldreich & Tremaine 1987). 
The mass-size relation of the GMCs is M c>c The frac- 
tal dimension of a volume fractal is D ~ 2.3 (Elmegreen 
& Falgarone 1996). The origin of this dimension could 
be turbulent diffusion in an incompressible fluid with a 
Kolmogorov velocity spectrum {D ^ 2 + for Ad oc 
see Meneveau & Sreenivasan 1990). 

The Kolmogorov theory applies for fully developed 
subsonic incompressible fluids. However, the ISM is su- 
personic and compressible. Only recently, 3D numerical 
studies of magneto-hydrodynamical and hydrodynami- 
cal turbulence in an isothermal, compressible, and self- 
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gravitating fluid indicated that the energy spectrum of 
supersonic compressible turbulence follows a Kolmogorov 
law (Mac Low 1999, Klessen et al. 2000, Mac Low & 
Ossenkopf 2000). The molecular clouds themselves are sta- 
bilized against gravitational collapse by the turbulent ve- 
locity field within them (see e.g. Larson 1981). 

Wada & Norman (1999) used high-resolution, 2D, hy- 
drodynamical simulations to investigate the evolution of a 
self-gravitating multiphase interstellar medium in a galac- 
tic disk. They found that a gravitationally and thermally 
unstable disk evolves towards a globally quasi-stationary 
state where the disk is characterized by clumpy and fila- 
mentary structures. The energy source of the turbulence 
in this system originate in the shear driven by galactic 
rotation and self-gravitational energy of the gas. The ef- 
fective Q parameter of the disk was found to have a value 
between 2 and 5. Without feedback the energy spectrum 
E(k) (X corresponds to a Kolmogorov law in two di- 
mensions, but changes into E{k) cx if stellar energy 
feedback is included (Wada & Norman 2001). This power 
law is expected if shocks dominate the system (Passot et 
al. 1988). Furthermore, Wada & Norman (2001) derived 
a driving length scale of ^200 pc for the model without 
stellar energy feedback. 

While the turbulent nature of the ISM is well estab- 
lished now, its origin and maintenance is still a matter 
of debate. This is of great importance, because turbu- 
lence can provide angular momentum transport in the disk 
of spiral galaxies. The evolution and the structure of an 
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accretion disk depends entirely on the effective viscosity 
caused by turbulence. 

Two ingredients arc necessary for a steady state tur- 
bulence: (i) a dynamical instability to generate and (ii) a 
steady energy input to maintain turbulence, 
(i) Five possible instabilities were put forward by different 
authors: 

— viscous instability (Lightman & Eardley 1974) 

— thermal instability (Pringle et al. 1973) 

— Parker instability (Parker 1966) 

— magneto-hydrodynamic instability (Balbus & Hawley 
1991) 

— gravitational instability (Toomre 1964). 

Because self-gravity is always present in the ISM, the grav- 
itational instability can operate in conjunction with all 
other forms of instabilities. As soon as the density rises 
due to such a combined instability and increases over the 
value of the ambient medium by a factor 2, self-gravitation 
should be the dominant force (Elmcgreen 1982). Moreover, 
the gravitational instability can produce a whole hierarchy 
of clumped structures inside the largest ones. Therefore, 
we will only discuss gravitational instability. 

(ii) Two main mechanisms are proposed to provide the 
energy input in order to maintain the observed turbulence 
of i^turb ~ 10 kms~^: 

— supernova explosions (see e.g. Spitzer 1968, McKee & 
Ostriker 1977) 

— galactic differential rotation (Fleck 1981). 

The supernova shocks, which can accelerate low mass 
clouds, are extremely ineffective in accelerating GMCs, 
because of the much larger mass to area ratio of GMCs. 
Jog & Ostriker (1988) applying directly the results of 
McKee & Ostriker (1977) found that supernova shocks 
can provide no more than 10% of the kinetic energy of 
the GMCs (however, Cui et al. (1996) suggested that un- 
certain background corrections could make this fraction 
larger). Therefore, we will only focus on energy input by 
galactic rotation. 

The effective viscosity VeS due to turbulence can be 
expressed in general as 



(1) 



where Wturb is the turbulent velocity on large scales, ^driv 
the driving wavelength, and ^ a constant. For continuous 
disks the most widely used form is the so-called a-ansatz 
introduced by Shakura (1972) and Shakura & Sunyaev 
(1973). They assumed that the turbulent velocity is lim- 
ited by the sound velocity (iiturb = Cg), the turbulent 
length scale equals the disk scale height Zariv = H, and 
^ — a < 1. Lin & Pringle (1987a) derived a prescription 
for the effective viscosity due to gravitational instability: 
Zdriv = H and ^ = where Q is the Toomre parameter 
(Toomre 1964). Recently, Duschl et al. (2000) suggested 
a prescription for the viscosity due to hydrodynamically 



driven turbulence at the critical effective Reynolds num- 

_ 1 

ber i?ecrit, where ^ = 1, i^turb = -Re^^f^ v^ot, and the driving 
wavelength is a well defined fraction of the local radius in 

the disk l^rw = I^^cTit ^■ 

For clumpy accretion disks, Goldreich & Tremaine 
(1978) and Stewart & Kaula (1980) elaborated models 
where the shear viscosity in a rotating disk is due to cloud- 
cloud interactions. Their prescription has the form: 



,, "^turb 

t^eS — 



1 



1 ^turb 



for r » 1, 



(2) 



where fl is the rotational angular velocity and r = f2 tint is 
the number of cloud interactions per rotation period with 
the interaction time scale tint between the clouds. Vertical 
hydrostatic equilibrium ~ I'rot/i'turb gives 



(3) 



where Vrot = is the rotation velocity. In this arti- 
cle we derive a prescription for the turbulent viscosity of 
comparable form with ^ = Re~^ in terms of the turbu- 
lent Reynolds number Re. All disk parameters are given 
as functions of Re, Q, and the mass accretion rate M. We 
apply our model to the Galaxy and derive a Galactic mass 
accretion rate. 

2. The basic picture 

We consider a gaseous accretion disk in a given gravita- 
tional potential $ which gives rise to an angular velocity 



n = ^JR-^'^■ The Toomre parameter (Toomre 1964) is 
treated as a free parameter 



Q 



Vtnrh K 



(4) 



with the restriction Q > 1, where G is the gravitational 
constant, S the average gas column density, and k the 
local epicyclic frequency. The disk consists of distinct self- 
gravitating clouds, which are orbiting in an external grav- 
itational field and have a velocity dispersion Av. These 
clouds might be embedded in a low density medium as 
long as Af/c™* < 1, where cf^* is the sound speed of the 
external medium. This is the case for the warm atomic 
or ionized phase (T ~ 8000 K) and the hot ionized phase 
(T ^ 10® K). The disk scale height is determined uniquely 
by the turbulent pressure Pturb = P'^^turb' where p is the 
averaged density in the disk. The clouds have a small 
volume filling factor <C 1. The transport of angular 
momentum is due to cloud-cloud interactions. These can 
be gravitational (elastic) encounters or direct (inelastic) 
collisions. Thus, gravitational scattering of the massive 
clouds off each other in a differentially rotating galac- 
tic disk gives rise to an effective "gravitational" viscos- 
ity. Jog & Ostriker (1988) determined analytically the 
one-dimensional velocity dispersion of such a system to 
be AuiD =5-7 kms~^ which is in good agreement with 
observations. We assume that energy is dissipated in the 
disk when clouds become self-gravitating. Turbulence is 
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generated by instabilities involving self-gravitation and is 
maintained by the energy input, which is provided by dif- 
ferential rotation via cloud-cloud interactions. Thus, the 
dissipative scale length Zdiss is equivalent to the size of 
the largest self-gravitating clouds. Since we assume clouds 
to form from local instabilities due to self-gravitation, 
which do not depend on the distance to the galaxy cen- 
ter, Q^const. Furthermore, we assume that the turbulent 
Reynolds number Re is also independent of the galactic 
radius. 

3. The equations 

3.1. The viscosity prescription 

In the inertial range of a turbulent medium kinetic energy 
is transferred from large scale structures to small scale 
structures practically without losing energy. So there is 
a constant energy flux from large scales to small scales 
where the energy is finally dissipated. Since the velocity 
dispersion (Aw ~ 10 kms~^) within the disk is more im- 
portant than the shear, there is no preferred transfer direc- 
tion. The turbulence is therefore assumed to be isotropic. 
In this case the similarity theory of Kolmogorov applies 
(see e.g. Landau & Lifschitz 1959). The assumption of a 
universal Kolmogorov equilibrium implies that the kinetic 
energy spectrum of the turbulence depends only on the 
energy dissipation rate per unit mass e and the character- 
istic size of the turbulent eddy I — -^j where k is the wave 
number. The kinetic energy E{k) is related to the mean 
kinetic energy in the following way: 

E{k)dk , (5) 



1 



2("(^)'> = 

where u is the velocity of the medium. Kolmogorov's the- 
ory yields 



E{k) = Ce3fc- 



(6) 



where C is a constant of order unity. Considering a 
schematic energy spectrum given by E(k) = 0, for k < 
/cdriv and for k > /cdiss and by Eq. |^ for fcdriv < k < fcdiss, 
one derives an expression for the dissipative length scale 
/diss — ^diss ^^'^ large scale velocity i^turb: 



/diss = (v^/^)^ , 

where v is the large scale viscosity due to turbulence; 



■^turb 



(7) 



(8) 



This leads to a relation between the two length scales and 
the turbulent Reynolds number, which is defined by Re = 



We will show in Sect. |^ that -Remacro ^ Re. In this work 
we will use the following viscosity prescription: 



1 

Re 



Wturb ^dr IV • 



(11) 



In our model the Reynolds number Re defined by Eq. ^ is a 
free parameter with the restriction Re > 1. This viscosity 
prescription will be compared to others in Sect. ^ 

It is formally equivalent to the expression for a clumpy 
accretion disk (Ozernoy et al. 1998) where the viscosity is 
due to cloud-cloud interactions (Eq. H). 

We interpret this equivalence as two different pic- 
tures for a turbulent self-gravitating medium. We pre- 
fer the point of view where the whole ISM (all phases) 
is taken as one turbulent gas which change phases (i) 
on turbulent time scales ttmb = ^turb/'f^turb and (ii) due 
to external processes (energy output by stars, i.e. super- 
novae, UV radiation by 0/B stars). Only near the mid- 
plane of the disk can the gas become molecular and self- 
gravitating which leads to a maximum GMC lifetime of 
At ^ Hmoi/t^turb ^ 10^ yr, where (i/moi ~ 100 pc Sanders 
et al. 1985). On the other hand the lifetime of a molec- 
ular cloud is approximately given by the crossing time 
At — dci/wturb, where dci is the cloud size. With a cloud 
size of 30 pc this results in At ~ 5 10^ yr. These values 
are comparable to the lifetimes of star-forming molecular 
clouds given by Blitz & Shu (1980). Since these lifetimes 
are about 10 times smaller than the expected mean col- 
lision time (Atcoii ~ 210^ yr Elmegreen (1987)), direct 
cloud-cloud collisions are very rare and are not important 
for the effective viscosity. 

Even with the limited lifetime of molecular clouds, 
gravitational interactions always take place between the 
continuously appearing and disappearing clouds. These in- 
teractions give rise to a collision term in the Boltzmann 
equation with an average collision time for gravitational 
encounters 



tG 



3^t*urb 



(12) 



(Braginskii 1965) where md is the cloud mass and rici is 
the spatial number density of the clouds. For molecular 
clouds the turbulent velocity is approximately indepen- 
dent of the cloud mass (iiturb oc to^j^ Larson (1981)), thus 



tG 



"-"'turb 



("^cl^-cl) 



with 



(13) 



Re'^ldiss ■ 



(9) 6.6 10"^ pc-3 / mlf dmd 



This turbulent Reynolds number is not equivalent to the 
macroscopic Reynolds number defined by e.g. Frank et al. 
(1992) 



Re 



(10) 



(14) 



where we have used the differential number density of 
molecular clouds given by Elmegreen (1987). Using inte- 
gration boundaries mi — 10** Mq, m2 = 10^ Mq and 
fturb^lO kms~^ one obtains tQ ^ 410^ yr. Assuming 
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a galactic radius R—10 kpc and a rotation velocity of 
t;rot=200 kms^^ gives t — Re ^ 80. 

Since the dissipative scale length is of the order of 
the largest self-gravitating clouds /diss ~ d^i ^ 30 pc, 
we can calculate the value of the driving wavelength 

3 

/driv = Re^ldiss — 800 pc. This is of the same order as the 
thickness of the distribution of the atomic gas (Kulkarni 
& Heiles 1987, Dickey 1993). We will show in Sect.| that 
in our model Idm ^ H. The scale height of the molecu- 
lar disk iJjnoi is smaller, 
molecular clouds. 



because of the limited lifetime of 



3.2. Non Kolomogorov energy spectra 

Wada & Norman (2001) found an energy spectrum of the 
form E(k) oc fc~^ in their 2D hydrodynamical simulations 
of a supersonic, compressible turbulence including stellar 
energy feedback. This power law corresponds to that of 
a shock dominated medium. In contrast to a compression 
dominated medium, the exponent of the power law for an 
incompressible medium is the same in two and three di- 
mensions (Passot et al. 1988). In this section we assume 
that E{k) oc in three dimensions In this case, an alter- 
native explanation for this exponent comes from the the- 
ory of intermittency (Frisch 1995). The open parameter 
of this model is the dimension D. In the case of intermit- 
tent turbulence the energy flux from scales ~ / to smaller 
scales is 



I 



, 3-D 



I H, 



= const 



driv 



(15) 



The local spatial average of the energy dissipation can 
then be expressed as 



I 



, D-3 



I 



^driv 



(16) 



cx fc ^ implies D — 2 and 
i.e. the local spatial average 



Within this framework E{k) 
consequently e; = e(|^^) , 
of the energy dissipation decreases with decreasing length 
scale I. The ratio between the driving and the dissipation 
length scale is then 



driv 

'diss 



= Re 



(17) 



instead of i?e4 in the case of Kolmogorov turbulence. For a 
Reynolds number Re ~ 50 the difference is about a factor 
2.5. Since the power law exponent of the ISM energy spec- 
trum is not yet well established, both possibilities must be 
considered. 

3.3. Angular momentum equation 

In a steady state accretion disk, the mass accretion rate 
is 



M = 27ri?S](-i;,ad) 



(18) 



where Wrad is the radial velocity. The angular momentum 
equation can be integrated giving 



i/E = -- 



M 



'dR' 



2ttR 

Furthermore, we use 



(19) 



(20) 



for the surface density of the disk, where p is the mass 
density at the midplane. 

3.4. Energy flux conservation 

The ISM turbulence in the disk is initiated by an insta- 
bility involving self-gravitation. The transport of angular 
momentum is due to cloud-cloud interactions giving rise 
to radial mass accretion. We assume that the necessary 
energy input to maintain the turbulence is the gravita- 
tional energy, which is gained when the ISM is accreted 
to smaller Galactic radii, i.e. the energy input is supplied 
by the Galactic differential rotation. Kinetic energy of the 
Galactic rotation is transfered to the turbulent cascade 
at the driving wavelength /driv and reaches the dissipative 
length scale /diss without loosing energy. The energy per 
unit time which is transferred by turbulence is 

E = p,. f V(T;Lb) df =^ pi^^A , (21) 

J Mriv 

where the integration is over the area / df = A. Thus, the 
energy flux per unit time and unit area is 



AE 

AtAA 



pu 



■^turb 



/dr 



(22) 



From standard fluid dynamics, the viscosity v generates 
dissipation in the disk at a rate 



1 ,v dn 

Mvroti^r- 

2tt dR 



(23) 



(see e.g. Pringle 1981), where Wrot ~ ft R. This is vahd 
for R ^ R^,, where {dO./ dR)ii=R, — 0. Conservation of 
energy flux leads thus to 



pv 



-^turb 



1 



/dr 



27r 



dR 



(24) 



and relates the mass accretion rate to the turbulence in 
the disk. 

3.5. The vertical pressure equilibrium 

Since we assume that the sound velocity of the ambient 
medium is smaller than the turbulent velocity dispersion 
of the clouds, the only pressure which counterbalances 
gravitation in the vertical direction is the turbulent pres- 
sure pturb — P^turb- ^c distinguish three cases for the 
gravitational force in the vertical z direction: 

1. Md{R) < 0.5 (H/R) M, where is the disk mass and 
M is the central mass (dominating central mass); 
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2. p < and Md{R) < Af,(i?), where p^/M^{R) is 
the stellar central density/disk mass within a radius R 
(dominating stellar disk mass); 

3. p > and 0.5 {H/R) M{R) < Md{R) < M(i?), where 
M{R) is the total mass enclosed within a radius R 
(self-gravitating gas disk in z direction). 

For these cases the gravitational pressure Pgrav has the 
following forms: 

1. pg,3^^ = pft'^H'^ (see e.g. Pringle 1981); 

2. Pgrav = ttGE*!], where E» is the stellar surface density 
of the disk (see e.g. Binney & Tremaine 1987); 

3. pg,a^ = ttGE^ (Paczyhski 1978). 

The equilibrium in the vertical direction is given by 



Pgrav — Pturb 



3.6. Global gravitational stability in z direction 



(25) 



The basic principles underlying the gravitational instabil- 
ity of a thin rotating disk can be found in Toomre (1964). 
A gaseous disk is locally stable to axisymmetric perturba- 
tions, if 



(26) 



where k = \ R 



4il2. Since in general Q. < k < 2^l, 



we will use the following equation: 

^ _ Wturb ^ 

^ ^ ttGE ■ 



(27) 



Multiplying the numerator and the denominator of the 
right hand side by R^ gives 



t'turb Mtot 
Wrot M^ 



(28) 



gas 



where Mtot/g^g is the total enclosed mass and the total 
enclosed gas mass at radius R. Thus for a given velocity 
dispersion is proportional to the ratio of gas mass to 
total mass. 



4. Results 

We have solved the set of equations Eq. |l^, Eq. |l^, Eq. 
Eq. H, Eq. ^ and Eq. for the three cases: 



1. dominating central mass, 

2. dominating stellar disk mass, 

3. self-gravitating gas disk in z direction. 

We will use 51' = dVl/dR with fi' < in the region of 
interest. 



4.1. Dominating central mass 

Pressure equilibrium (Eq. p5| ) leads to 



(29) 



which is equivalent to the hydrostatic equilibrium where 
the sound speed is replaced by the turbulent velocity dis- 
persion. 



From energy flux conservation (Eq. 24) and the angu- 
lar momentum equation (Eq. |l9|) it follows 



H 



(-R-^) 'driv 



(30) 



Thus, the driving wavelength is proportional to the disk 
height /driv oc H and the viscosity prescription reads v ^ 
Re'^vtuihH, which is formally equivalent to the viscosity 
prescription for clumpy disks of Goldreich & Tremaine 
(1978), Stewart & Kaula (1980) (see Eq. |), and the a- 
viscosity (Shakura & Sunyaev 1973) with Wturb = Cs and 
Re-^ = a. 

Inserting Eq. in Eq. and using p = T./H gives 



ttQ G 



(31) 



If we define the critical cloud density for stability against 
tidal shear pcrit = 7r~^ri^G~^, the Toomre parameter Q 
represents the ratio between the critical density and the 
volume averaged density in the disk Q = Pcrit/p- 

Inserting the viscosity prescription (Eq. |ll|) into the 
angular momentum equation (Eq. |l9| ) leads to 



H 



{^GReQMR{-n')n-^) 



With S = p/H we obtain 



pH 



1 



237r 

Inserting Eq. pSl into Eq. [19 



G~^Re^Q-^M^R^{~n')^ni 



gives 



1 

2l 



GiRe-^Q^M3R-3(-n')-3n3 



(32) 



(33) 



(34) 



The turbulent velocity can be calculated using the pres- 
sure equilibrium (Eq. p5|): 



«turb 



.1 



GReQ M R{-n')n- 



(35) 



If il' ~ —il/R, the viscosity is inversely proportional to 
the angular velocity v ex ^l^^ and the turbulent veloc- 
ity is constant. In the case of a Keplerian velocity field 
due to a point mass f cx R^^^. The behaviour of the disk 
parameters (vturb, 'driv, S, H, p, and ly) are shown for a 
constant rotation curve (solid line) and a rising rotation 
curve Vrot oc VR (dashed line) in Fig. [|. The free param- 
eters (5 = 1, Re=50, and M=10~^ Moyr~^ were chosen 
such that the disk parameters fit the observed values for 
the Galaxy (see Sect. 0). 
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Fig. 1. Disk parameters for the case of a dominating cen- 
tral mass. Solid lines: constant rotation curve. Dashed line: 
rising rotation curve frot oc y/R. If Q=l the behaviour of 
these parameters is the same in the case of a selfgravitat- 
ing gas disk in z direction. 



4.2. Dominating stellar disk 

Pressure equilibrium (Eq. |5|) together with f^^^ 
TzGT... R leads to 



(36) 



Inserting Eq. |3^ in the equation for the Toomre param- 
eter (Eq|2^) and using the viscosity prescription (Eq. |ll|), 
the angular momentum equation (Eq. |l^) and energy flux 
conservation (Eq. 24) gives 



(37) 



In the case of 57' ~ fi/i?, the turbulent velocity is inversely 
proportional to the square root of the rotation velocity 

-1/2 

Inserting Eq. ISTI into Eq. 



1 



H = ^GReQMR-^{-n')n-^ 



gives the disk scale height 
(38) 



Eq. |ll|, Eq. and Eq. H together with Eq. |3| lead to 
the following expression for the driving wavelength 



^driv — (^^') ^ ^ ^ 



(39) 



From the same equations the density, column density, and 
viscosity read 



/2 

TT 



(40) 



Fig. 2. Disk parameters for the case of a dominating stel- 
lar disk. Solid lines: constant rotation curve. Dashed line: 
rising rotation curve v^ot oc y/R. 



1 



V27r 



G"5Q-2i?e2M2(_rj')! 



and 



= V2G^Q^ReriMiR'^{-n')-^n 



(41) 



(42) 



If Q' ~ — fl/R, the driving wavelength is proportional 
to the galactic radius Z^riv R. Therefore, in the case of a 
dominating stellar disk, one obtains for the viscosity pre- 
scription (Eq. |ll|) I' ^ Re^^ wturb R- In this case, the vis- 
cosity is proportional to the galactic radius divided by the 

— 1/2 

square of the rotation velocity ly cx RVj.^^ . Furthermore, 
if the rotation curve is flat, i.e. Urot=const., u oc R. This 
behaviour can be observed in Fig. || where the disk pa- 
rameters (wturb, ^driv, S, H, p, and v) are shown for a 
constant rotation curve (solid line) and a rising rotation 
curve Urot oc Vi? (dashed line) for (3=1000, i?e=50, and 
A/=10^^ Moyr^^. Re, and M are the same as for the 
Galaxy (see Sect. |^) and the value of Q is chosen such 
that the turbulent velocity dispersion corresponds to that 
measured of an SO galaxy (~50 kms~^ D'Onofrio et al. 
1999). 

4.3. Self-gravitating gas in z direction 

The pressure equilibrium (Eq. p5|) together with the equa- 
tion for the Toomre parameter (Eq. leads to 



P - 



1 

7rg2 G 



(43) 



The square of the Toomre parameter thus represents the 
fraction between the critical density with respect to tidal 



B. VoUmer and T. Beckert; Turbulent viscosity in clumpy accretion disks 



7 



shear and the volume averaged density in the disk = 
Pcrit/p- We want to recall here that for a dominating point 
mass Q = Pcrit/p- 

Inserting Eq. ^ into the pressure equilibrium equation 
(Eq. GH) gives 



H = QVL ^ Wturb 



(44) 



Then, using the angular momentum equation (Eq. |l9|) , 
it follows 



(45) 



From the energy flux conservation (Eq. |2J) one obtains 
A 



"turb ^ [-GQ^ ReM R {-Q') n 



-l^ 3 



(46) 



For n' ~ ~^l/R, the turbulent velocity is constant. This 
is what is observed in the nearby galaxy NGC 3938 (van 
der Kruit & Shostak 1982). 

Inserting this expression in Eq. nB leads to 



tdriv 



(47) 



For fl' ~ —Q/R, H/l^nv — ■ Thus, the viscosity pre- 
scription yields v — Re~^Q~'^VtuihH . Here, Re is not the 
macroscopic Reynolds number Remacm = VmtR/v (see 
Sect. In the case of i?e = 1 this is the prescription 
given by Lin & Pringle (1987a). 

Inserting Eq. ^ into Eq. ^ leads to the disk height 



1 



H = ^QiRe^G^M^R^{-n')^n- 

23 



(48) 



With S = p/H one obtains the following expression for 
the column density: 



E = -^Q-^^Re^G-iM^R^{-n')^n^ 
23 7r 



Eq. |Tl|, Eq. |46[ Eq. and Eq. || lead to 

V = Xo^Re-^G^MiR^U-n'y^n^ . 
2§ 



(49) 



(50) 



In the case of Q,' ~ —Q./R, the column density is pro- 
portional to and the viscosity is inversely proportional to 
the angular velocity: E cx 51, cx fl^^. If in addition the 
rotation velocity is constant, then v (x R. 

Comparing this result with Eq. ^ shows that both 
models with a gravitational potential $ due to an ex- 
tended mass distribution give the same viscosity prescrip- 
tion if the rotation velocity Urot = ^i?d$/di? is constant. 

On the other hand, the dominating central mass model 
and the self-gravitating gas disk model have the same an- 
alytical solutions for Q = 1 (see Fig. |l|). 



5. Discussion of the results 

5.1. Comparison with previous models of turbulent, 
selfgravitating gas disks 

Models of turbulent, self-gravitating gas disks can be di- 
vided into two different approaches: (i) The disk is as- 
sumed to be quasi continuous and at the edge of fragmen- 
tation {Q ~ 1); (ii) The disk is already clumpy and the 
viscosity is due to cloud-cloud interactions {Q > 1). 

We will first discuss the quasi continuous approach. 
Paczyhski (1978) investigated a self-gravitating disk with 
a polytropic index of 7 = | which corresponds to a radi- 
ation pressure dominated disk. Furthermore, he assumed 
the disk luminosity to be at the Eddington limit. With 
the acceleration due to a central mass M and the disk 
gas surface density 'S: g — GMR^^z + 2ttGYi, he ob- 
tained S cx R^^ and thus v cx R^ for a constant ro- 
tation curve. Lin & Pringle (1987a) proposed a viscos- 
ity prescription based on the Toomre instability criterion 
(Toomre 1964): v = Q-'^H'^R. They showed that under 
certain conditions this prescription allows a similarity so- 
lution E cx R^ 2 . Their viscosity prescription can be gen- 
eralized v cx E°i?-'' with a, /? > or V fx ^ for 
self-gravitating disk in z and R direction (Saio & Yoshii 
1990). Shlosman & Begelman (1989) also considered a disk 
at the edge of selfgravitation, i.e. Q ^ I. They obtained 
E cx t'turb'^'rot^"^ oc R~^ for coustaut turbulent and rota- 
tion velocities. On the other hand, Duschl et al. (2000) 
made a completely different ansatz: v cx Pv^otR, with 
/3 <C 1. They found E cx R~^ and ;/ cx i? for a disk with 
Vrot = const. 

On the other hand, using a clumpy disk model Silk & 
Norman (1981) derived a viscosity prescription in assum- 
ing that the cooling time for cloud-cloud collisions tcooi = 
^o(^^turb'7)~^ equals the viscous time scale = R^v~^ for 
all radii R, where 77 is the fraction of cloud kinetic en- 
ergy radiated in a collision and is the cloud mean free 
path. This leads to v = vo{R/Rq), where vq — y^-Ro'f^turb 
and Rq is a characteristic length scale. Lynden Bell & 
Pringle (1974) showed that for this viscosity prescrip- 
tion E cx (i?t)~^ exp(— ^i?/t), where ^ is a constant and 
t is time. Later Shlosman & Begelman (1987) also used 
^cou ^ tvisc- Ozernoy (1998) and Kumar (1999) suggested 
a viscosity prescription based on the coUisional Boltzmann 

Tv'^^^^n, where r = (tcoii^^)"^- 



equation (see Sect. 3.1) v 



In this section we will only discuss the self-gravitating 
gas disk in z direction, because only for this set of equation 
Q ~ 1, i.e. the gravitational instability is active. This point 
will be further discussed in Sect. 



5.2 



Whereas the above cited viscosity prescriptions (ex- 
cept Duschl et al. 2000) use Zdriv = H, the driving wave- 
length Zdriv is a priori a free parameter in our set of equa- 
tions. It is then defined by the energy flux conservation 
(Eq. I2J). In terms of the standard accretion disk equa- 
tions (see e.g. Pringle 1981), the "thermostat" equation, 
i.e. the viscous heating is radiated by the disk , is replaced 
by the turbulent energy flux equation (Eq. |2j) , where the 
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energy gained through accretion and differential rotation 
is transported by turbulence to smaller scales where it is 
dissipated. The addition of this energy flux conservation to 
the "traditional" set of disk equations leads to Zdriv = H . 
Once this relation is established, the viscosity prescription 
has the form of that derived from the collisional Boltzman 
equation with Re — tcow^ (see Sect. 3.1). On the other 
hand, our viscosity is a factor Re~^ smaller than that of 
Lin & Pringle (1987a) for Q - 1. 

Thus, a big difference to previous models is that we 
do not use a disk "thermostat" , i.e. that the viscous heat- 
ing is radiated by the disk. Radiation is not included in 
our model. Duschl et al. (2000) noticed that the a viscos- 
ity together with a pressure term due to self-gravitation 
leads to a constant sound velocity. They declared this un- 
physical within the framework of a disk "thermostat" . On 
the other hand, if one replaces the sound velocity by the 
turbulent velocity this is exactly what is observed in spiral 
galaxies. Since t^turb = const we have H/R — const., and 
consequently v oc R and S oc R~^ . 

5.2. Why can the Toomre parameter be greater than 
1? 

Whereas the two approaches of a continuous and a clumpy 
disk can be unified, this can not be done for a disk with 
Q > 1. In this case the disk is not globally gravitational 
unstable (Toomre 1964). The turbulent velocity is so high, 
or the surface density is so low that the disk is globally 
gravitationally stable. However, we can imagine two pos- 
sibilities for the formation of Q > 1 disks: 

(i) In Sect. |^ we will show that the mass accretion rate 
within the disk is much smaller than the star formation 
rate. Thus, starting from a Q ~ 1 disk with star formation, 
Q will increase with time (if there is no external mass 
infall). This could be the case for SO galaxies. 

(it) The gas which falls into the given gravitational 
potential is already clumpy. This could be the case for 
the Circumnuclear Disk in the Galactic Center (VoUmer 
& Duschl 2000). 

For these scenarios the above models for a dominating 
central mass and a dominating stellar disk are valid. We 
are the first to solve these sets of equations. 

We will now show that Q = const > 1 can be un- 
derstood in terms of star formation. With the critical 
density for tidal disruption pcrit = ^'^/{ttG) and using 

= Re~'^T,n (Eq. |l|) one obtains: 



P* = 



1 Pcrit 

Re Q in 



where tu ~ H/vtuvh is the vertical crossing time. Thus, 
Q = const > 1 together with Re — const > 1 means 
that the star formation rate is proportional to the critical 
density with respect to tidal shear divided by the vertical 
crossing time. 

The detailed comparison between observations and our 
model will show if Q > 1 models are useful to describe 



clumpy disks with a low gas mass or high velocity disper- 
sion. 

5.3. Energy dissipation due to selfgravitation 

If turbulence is dissipated by self-gravitation, the energy 
dissipation rate due to self-gravitation of the disk in z 
direction egg must equal the constant, turbulent energy 
dissipation rate per unit mass eturb- 

The energy due to self-gravitation of a gaseous disk is 



AE ^pAAH ^ ttGY? AAH = irGT^lVL^^H , 



(52) 



where is the disk surface. The growth rate of the gravi- 
tational instability in z direction is approximately the free 
faU time At - ^JH/{GT.) = ^{Gp)-^. Thus, one obtains 
for the energy dissipation rate per unit mass due to self- 
gravitation 



AE 



-kGYjH^Gp . 



AM At 

Using Eq. ^ and Eq. ^ it can be shown that 



V 



Cturb 



V. 



turb 



turb 

H 



(53) 



(54) 



In order to estimate the mass of a self-gravitating cloud, 
we use AE ^ N^^M^^G /Idiss, where A'ci is the number of 
selfgravitating clouds and Mc\ is the cloud mass. With the 

y(Gppi$v this leads 



free fall time of the clouds At 

to 



^ 3 _ 
AM At 5 



Gp$v = etu 



(55) 



where is the volume filling factor (see Sect. p.l\) and 
AM = N ciMci. Using the expression for $v derived in 
Sect. |t| and Eq. || gives 



Mr, 



-G-'Re-^Q-'n 



^turb 



(56) 



210-8 yr-i, and 
210^ Mf.. This is 



With Q = 1, Re=50, ft = 
t^turb=10 kms"^ one obtains Md 
consistent with the mass of GMCs. 

In the present paper we have assumed that the poten- 
tial energy gained through accretion is dissipated locally. 
This is true only in the absence of radial energy transport, 
which can be checked a posteriori. The energy dissipation 
rate due to self-gravitation of the disk in z direction is 



(51) 



AE t:GY?H 



AAAt tn 
The radial energy flux is given by 



(57) 



(58) 



where — —M/{2ttRY,) is the radial velocity and e 
the specific internal energy. For supersonic turbulence in 
clumpy disks e is dominated by the specific kinetic energy 
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v'^^j.y^. Since iiturb = const, de/dR ^ 0. With v-^ = const 
one obtains 

,5, - -2E<.,^;,i?-i = -2pvl^^^v,HR-^ . (59) 
Using <H = H/vt^i\y and Eqs. ^ one obtains 



5r 



2^3TrQ3Re-3G3Msv. 



2 

rot 



(60) 



With Q=l, i?e=50, and M = IQ-^ Mgyr"! (see Sect, 
this leads to 



^sg (wrot[kms ^])2 



< 1 



(61) 



The radial advective energy flux is much smaller than the 
energy gain by mass accretion and can therefore be ne- 
glected as long as Eq. ^ holds. 

We can conclude here that the energy flux transported 
to small scales by turbulence has the same order of mag- 
nitude as the energy dissipation rate due to selfgravita- 
tion of the gaseous disk. Moreover, the potential energy 
gained through accretion is dissipated locally at small 
scales. Since our model flts observations for a reasonable 
choice of parameters, we conclude that, based on energy 
flux conservation, it is possible that turbulence is gen- 
erated and dissipated through gravitational instabilities, 
and maintained by the energy input due to mass inflow 
and differential rotation (Sect. [3.4[ ). 

6. Implications 

6.1. Volume filling factors 

In this Section we compare the crossing time of a turbulent 
cloud to the gravitational free fall time in order to derive 
an expression for the volume filling factor as a function 
of the Reynolds number Re and the Toomre parameter 
Q. We assume turbulence with a Kolmogorov spectrum 
for Zdiss < I < 'driv This implies: /diss ~ Re~^ldriv and 
Vci = Re~3VtuTin where Vd is the turbulent velocity of a 
cloud of size /diss- 

The characteristic turbulent time scale of clouds whose 
size is comparable to the dissipative length scale is 



tl = /diss/Wcl ■ 

The gravitational free fall time is given by 



iff = 



37r 



(62) 



(63) 



where pd is the density of a single cloud, which is related 
to the overall disk density p by the volume filling factor 

(l>v- pel = (l>Y^P- 

The clouds are self gravitating for t\ = iff. Inserting 
the expressions given in Sect. leads to: 



1. dominating central mass: 

^Q-^Re-^R-^{^n') 



Dy 



dominating stellar disk mass: 

= ^G-iRe-iQ-'iiriR-\^n')-"^n^'' 

= ^Rl~^Q-^R-\-^')-'n\^f. 
self-gravitating gas disk in z direction: 
Q-^Re-^R-^{-Vl')~^n^. 



32 



For a non Kolmogorov spectrum E{k) (x k ^, the 
volume filling factors have to be multiplied by Re~^: 



6.2. Bondi-Hoyle accretion limit 

We will now consider the limit (j)y = Q = Re ^ 1, which 
means A/gas ~ Mtot and wturb ~ Wrot- This implies H ~ 
R, i.e. this is the limit for a spherical configuration. It is 
assumed that dfl/dR ~ —Q./R. The viscosity prescription 



(Eq. 
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(Eq. 24) gives 



together with energy fiux conservation equation 



1 3 _ Af n 



(64) 



Setting Wturb = Vrot, Re = 1 and using v'^^^ = MG/R leads 
to 



M = 



2ttG^M^P 



(65) 



This is equivalent to the Bondi-Hoyle accretion rate 
(Bondi & Hoyle 1944), which is given by 



M, 



B-H 



2.5ttG^A'Pp 



(66) 



turb 



for Vturb = Vr. 



6.3. The molecular fraction 

In order to derive an expression for the molecular fraction 
of the gas in the disk, we compare the crossing time of the 
turbulent layer tturb and the H-H2 transition time scale 
tii2 = a/p (HoUenbach & Tielens 1997). We define the 
molecular fraction here as fmoi — ^turb /%2 • 

For the three different cases we obtain using a 
Kolmogorov spectrum: 

1. dominating central mass: 

/,noi = l/i7ra)Q-^Re-iG-^R-^{-n')-'^^ 

2. dominating stellar disk mass: 

/„,oi = 2/{7ra)G~^Q"^Re~iM-\-n')-^n^ = 

= 2/(™) G-^Q-^Re-lM-^R-^i-n'y^vf,^n\ 

3. self-gravitating gas disk in z direction: 

= l/{7Ta)Q-^Re-^G-^R-H-n')-^n^ 

For a non Kolmogorov spectrum E{k) cx k^^, the 
molecular fractions have to be multiplied by Re^^: 

fnoiiK. J^f^~h fl^- 

J-mo\ — - J-moV 

In the following we will only treat the case of a self- 
gravitating gas disk in z direction, because it represents 
the most realistic description of normal fiel d sp iral galax- 
ies. We will motivate this choice in Section 6.4. 
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For a galaxy with a constant rotation velocity the 
molecular gas surface density is given by 

= — Q-^Re-iG~Hlin^ . (67) 

The total integrated molecular gas mass is then 

= TrSn^i?^ = -^—Q-^Re-iC'^Mh'^^^ , (68) 

which does not depend on the galactic radius. 

With these results we can now calculate the ratio of 
molecular to atomic gas in spiral galaxies 

Mh, /AfHi = Mh, / ( Af - Mh J , (69) 
where Af^"^ = TrSi?^ is the total gas mass. 

6.4. Star formation and the Tully-Fisher relation 

The star formation rate per unit area is usually de- 
scribed by a Schmidt law (Schmidt 1959) of the form 

(X Ef,, , (70) 

where ^ 1.4 (Kcnnicutt 1998). An alternative descrip- 
tion of the star formation rate was proposed by Elmegreen 
(1997) and Silk (1997). They suggested that it scales with 
the ratio of the gas density to the average orbital time 
scale 

±, = 0.017Egasf^ . (71) 

This parameterization provides a fit that is nearly as good 
as the Schmidt law. 

On the basis of our model, we suggest that the star 
formation time scale is tsf = Reil~^, which is the in- 
teraction time scale between the clouds. This corresponds 
to the time scale of gravitational encounters between the 
clouds. We thus find Re ^ 60, which is only a factor 2 
higher than our estimate based on the disk height and the 
cloud size. In the case of a constant rotation velocity our 
star formation prescription yields in the case of selfgravi- 
tation in z direction: 

E* = -^Q-i Re-^G-i Klin"^ . (72) 

23 TT 

Thus both, the star formation rate per unit area and the 
molecular gas density, are proportional to the square of the 
angular velocity oc fi^ cx fi^. Both profiles have 
thus the same dependence on the galactic radius. This is 
what is found when comparing CO and Ha, radio contin- 
uum, B band luminosity profiles of most of the Virgo clus- 
ter galaxies observed by Kenney & Young (1989). More 
recently, Rownd & Young (1999) combined Ha imaging 
and CO-line observations to derive radial star formation 
efficiencies for a sample of cluster and field spiral galax- 
ies. They found that E^,/Eh2 = const, which is consistent 
with our findings. It is worth noting, that in the frame- 
work of our model this proportionality only exists between 
the star formation rate and the molecular gas mass and 
not the total gas mass. 



Whereas the star formation laws of Kennicutt (1998) 
or Silk (1997) depend explicitly on the gas surface density, 
our prescription only depends on the angular velocity f2. 
The proportionality factor depends on M and the free pa- 
rameters Re, Q. The results of Rownd & Young (1999) im- 
ply that the combination of these parameters or all param- 
eters individually are the same for all spiral galaxies. The 
prescription of Silk (1997) is equivalent to ours E* = ^Ef2. 
We identify ^ = i?e. The observationally derived ^ is con- 
sistent with the definition of Re = (Zdiiv/'diss)' • For the 
comparison between the Kcnnicutt and Silk law we refer 
to Kennicutt(1998). 

The integrated star formation rate is 

= ^{^)'m = ^Q-iRe-iG-lMhl, . (73) 

2 Vtnrh 23 

Thus the integrated star formation rate is proportional to 
the square of the rotation velocity as the total amount 
of molecular gas. Observationally, Devereux & Hameed 
(1997) found that the ratio of star formation rate to molec- 
ular gas (Lfir/A^(H2)) is independent of the morpholog- 
ical type of a spiral galaxy, thus Af * oc AfH2 ■ 

The blue Tully-Fisher relation is dominated by light 
associated with current star formation. In general, the 
galaxy luminosity is proportional to the rotation veloc- 
ity at a certain power L oc ii^ot- Compilations of Tully- 
Fisher data for available samples find that 7=2.1-2.2 for 
the B band with a systematic increase with increasing 
wavelength (Silk 1997). Our model would ideally predict 
Af» oc w^Q^., i.e. 7 — 2. Since the B band luminosity is pro- 
vided by a mixture of stars of different ages, we would ex- 
pect that 7 is slightly higher than predicted by our model. 

We will now calculate the integrated star forma- 
tion rate for the case of a galaxy with dominating stel- 
lar disk. Inserting the expression for the volume fill- 
ing factor for a constant rotation velocity disk (/)v — 
^ Q~^Re~^{viot/vtuih)^ and Eq. ^ into the expression 
for the integrated star formation rate A/» = TTl]R'^Re~^n 
yields 

Af, = ^^Af = (^i?eQ$v)3M . (74) 

t^turb 4 

Taking (/.y 10"^ Re = 50, M = 10"^ Mgyr"!, 
and Urot = 200 kms^^ gives utuib 50 kms~-^ and 
Af* ~ 410"^ Mqyi'^. Thus, the integrated star for- 
mation rate is comparable to the mass accretion rate. 
Such low integrated star formation rates are only ob- 
served for SO galaxies (Kennicutt 1998). The Toomre pa- 
rameter for such a disk is Q ~ 710'^. For ipy — 1, 
Q Re^ — ^{vmt/vtmhY ^ 10^- This is not the case for nor- 
mal field spirals. Thus, for most of the field spiral galaxies, 
the case of a self-gravitating disk in z direction applies. 

6.5. Low surface brightness galaxies 

Low surface brightness (LSB) galaxies have low den- 
sity/column density gas but a high gas fraction A/gas/Aitot 
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(de Blok 1999). These LSB galaxies can be either dwarf 
galaxies with small rotation velocities (see e.g. Walter & 
Brinks 1999) or large galaxies with high rotation velocities 
(see e.g. Pickering et al. 1997). Both have small angular 



Q= -.00000 Re= 50,0000 



velocities 51 = v^-ot/R- Since om' model yields S, 



gas 



and p cx il^, a small angular velocity implies a low gas col- 
umn density. Thus, the observed small gas density /column 
density of LSB galaxies can be partly due to their small 
angular velocity. 

7. Application to the Galaxy 

As motivated in Sect. |6.4| the model of a vertically self- 
gravitating gas disk describes the gas disks of spiral galax- 
ies. We will use the local parameter of the ISM in the so- 
lar neighbourhood given by Binney & Tremaine (1987): 
p = 0.042 Mqpc-3, E = 5.3 Mgpc-^, = 2.62 10^^ yr'^ 
and i?=10 kpc. Furthermore, we have shown in Section |3.l| 
and Section 6.4 that Re ^ 60-80. The rotation curve is 
assumed to be constant. Eq. 43 leads to Q ~ 1. On the 
other hand, the definition of the Q parameter (Eq. 
would yield Wtmb ~ 3 kms~^. This equation would thus 
require Q ^ 3 in order to match the gas velocity disper- 
sion of Wturb=10 kms^^ found by van der Kruit & Shostak 
(1982). We will here adopt 

Q = i, 

Re = 50, 

a = 10*" yrMgpc"^ (this corresponds to a molecular tran- 
sition time scale of tmoi ~ 5 10^/n yr, where n is the par- 
ticle density), 
i?=10 kpc, and 
t^tuib=10 kms^^ 

Eq. 1^ then leads to M 10"^ Mgyr-^ Thus, we 
suggest that the Galactic mass accretion rate is M ~ 
10"^ Mgyr"^. Fig. ^ shows the derived turbulent velocity, 
driving wavelength, gas column density, disk scale height, 
gas density, viscosity, molecular fraction, and volume fill- 
ing factor for three different rotation velocities. For the 
model of the Galaxy, we have adopted 
i'rot=250 kms~^. 

This results in a total gas mass of Afg°* — 5.4 10^ M©, 
a total molecular gas mass A/hs = 1.710^ Mg, a total 
atomic gas mass Afni = 3.710^ M©, and M^i^l^m = 
0.3. ScoviUe & Sanders (1987) give Mm ~ 2.5 10^ Mq, 
A/hj ~ 0.7 10^ M0 (we have used the conversion factor 
X = 10^0 cm-2(Kkms"i)-i of Digelet al. 1996) and thus 
M-a^/M-Hi ~ 0.3. Our model calculations are consistent 
within a factor of 2 with the absolute values of the total 
gas mass and the observed ratio of molecular to atomic 
gas in the Galaxy. Figure ^ shows the radial profiles of 
the total gas (dotted line), the molecular (solid line), and 
the atomic (dashed line) gas surface density for Uj-ot = 
250 kms""'^ and a — 710^ yrMopc^. These profiles well 
resemble that found by Kenney & Young (1989) (their 
Fig. 12). Honma et al. (1995) examined the Hi and H2 
radial distributions of four face-on spiral galaxies. They 
found that the gas phase transition occurs suddenly within 
a narrow range of radius. This transition is explained using 
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Fig. 3. Radial profiles of the rotation velocity, angular ve- 
locity, turbulent velocity, driving wavelength, gas column 
density, disk scale height, gas density, viscosity, molecular 
fraction, and volume filling factor for three different rota- 
tion velocities (solid line: v^t =250 km/,s~"'^, dashed line: 
^'rot =175 km/,s^^, dotted line: Wrot =100 km/,s~^). 



a theory whose main parameters are the ISM pressure, the 
UV radiation field, and the metallicity. It was found that 
the metallicity gradient is most crucial for the formation of 
a narrow range transition. Our model (Fig. ^ also shows 
a sharp phase transition. For a given i?e, Q, and Af, the 
location of this transition region depends mainly on the 
formation time scale a of H2 from Hi. If a depends on 
metallicity (and that is what one would expect), the two 
models are in good agreement, because radiation, which 
is not included in our model, does not play a major role. 

The disk height equals the driving wavelength H ~ 
^driv ~ 400 pc, and the viscosity is ^ 8 10"^ pc^yr~^ = 
2 lO^'^cm^s^^. Thus, the macroscopic Reynolds number is 
-Remacro = v^otR/v ~ 3 10^ ^ Re. This value of the driv- 
ing wavelength is only a factor of 2 larger than that found 
by Wada & Norman (2001) in their 2D hydrodynamical 
simulations. 
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Fig. 4. Radial profiles of the total gas (dotted line), the 
molecular (solid line), and the atomic (dashed line) gas 
surface density. 



The derived mass accretion rate is much smaller than 

" -1 



the star formation rate (Eq. 73 yields M* ~ 3 M0yr 



compared to ~ 3 — 5 Moyr"-^ as suggested by the 
observations of other Sbc spirals; Tammann et al. 1994; 
Cappellaro et al. 1997). The current viscous time scale 
~ R^/i/ ^ 10^^ yr is larger than a Hubble time. Since 
the viscous time scale is given by 



i?2 



(75) 



viscous evolution has taken place at a time when the mass 
accretion rate was higher and/or the disk was smaller 
and/or rotation velocity was smaller. After this evolution 
the gas distribution reached an equilibrium state that is 
described by our model. Lin & Pringle (1987b) and later 
Saio & Yoshii (1990) have shown that one obtains an ex- 
ponential stellar disk if ti, ~ t*. In our model this is only 
the case for a dominating stellar disk, where M« ~ M. 
However, we have to stress here that M is the mass ac- 
cretion rate in the disk. If there is also external accre- 
tion present, this formally increases i*, which can lead to 



8. Morphological type versus mass 

In this Section we investigate the link between the mor- 
phological classification of galaxies and their mass on the 
basis of our model for vertically selfgravitating gas disks. 
In general, disk galaxies of later types (Sc-Scd) are smaller 
and less massive than galaxies of earlier type (Sa-Sab) 
(Roberts & Haynes 1994). Some of the properties which 
lead to the morphological classification could thus be due 
to the galaxies' gravitational potential as proposed by 
Gavazzi et al. (1996). 
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Fig. 5. Ratio of molecular to atomic gas mass as a func- 
tion of the galaxy rotation velocity. 

8.1. The velocity dispersion and the fraction of gas to 
total mass 

It is a surprising fact that galactic gas disks all have disper- 
sion velocities between 8 and 10 kms^^ (see e.g. Freeman 
1999). Assuming a constant rotation velocity Eq. ^ im- 
plies that Re M is approximately constant for spiral 
galaxies of all types and all luminosities. Furthermore, if 
the galaxy forms still actively stars Q should be r^l. Thus 
Re M can only vary by a factor ^10. 

From Eq. || we find Mgas/Mtot = Q'^vturh/vrot- 
Since the gas velocity dispersion is almost constant, 
Mgas/Mtot oc Galaxies with high rotation velocities 
have thus smaller gas to total mass ratios. This trend is 
also observed in the Hubble sequence: early type spirals 
have larger rotation velocities and smaller gas to total 
mass ratios (Roberts & Haynes 1994). 

8.2. The fraction of molecular to atomic gas 

Young & Knezek (1989) have shown that the ratio of 
molecular to atomic gas mass in spiral galaxies increases 
from early type galaxies to late type galaxies by a factor 
~ 10. We have calculated the ratio 



Ah 



H2 



/mol S R'^ 



/mol 



MhI TT S i?2 „ /mol S i?2 1 - /„iol 



(76) 



We used the following parameters: Q = 1, Re — 50, 
a — 10^ yrMopc"^, and i?=10 kpc. The result can be seen 
in Fig. There is almost a factor 10 between the molec- 
ular to atomic gas ratio of a galaxy with v^-ot = 50 kms~^ 
and a galaxy with Vrot = 300 kms""'^. Of course, galaxies 
with small rotation velocities tend to be smaller, but their 
Hi diameter should be larger, justifying i?=const. as a first 
approach. Furthermore, we expect that there is an influ- 
ence of the molecule formation time scale a on the ratio 
between molecular and atomic gas mass. This time scale 
should depend on the galaxy's metallicity. Thus the ob- 
served trend might be a mixture of (i) higher metallicities 



B. VoUmer and T. Beckert: Turbulent viscosity in clumpy accretion disks 



13 



(and therefore smaller a) of early type spiral galaxies and 
(ii) higher rotation velocities of early type spiral galaxies 
compared to late type spirals, i.e. that early type galaxies 
are more massive. 

8.3. Star formation rate per unit mass 

The integrated star formation rate times the rotation pe- 
riod divided by the total mass should be comparable to the 
ratio of FIR luminosity to H band luminosity LpiR/ifj. 
Using Eq. ^ and w^qj ^ MtotG/R, we obtain for our 
model 



This means that rapidly rotating galaxies have smaller 
Lpm/Ln- Devereux & Hameed (1997) found indeed that 
galaxies of later type have smaller Lfir/^h- Since the 
FIR luminosity is assumed to be a measure of recent star 
formation and the H band luminosity a measure of the to- 
tal mass, (Af*f2~-'^)/Aftot decreases with increasing mass. 
Moreover, Gavazzi et al. (1996) found a trend of smaller 
Ha equivalent width for more luminous (H band) galax- 
ies. This also translates into a decreasing (M,r2~^)/Aftot 
with increasing galaxy mass. Thus, again, a part of this 
observed trend might be due to the fact that early spirals 
have higher rotation velocities than late type spirals. 

9. Summary 

We have analytically investigated the equilibrium state 
of a turbulent clumpy gas disk. The disk consists of dis- 
tinct self-gravitating clouds, which are be embedded in a 
low density medium, and evolve in the fixed gravitational 
potential of the galaxy. The disk scale height results from 
the balance between the gravitational acceleration and the 
pressure due to the turbulent velocity dispersion of the 
clouds. Gravitational cloud-cloud interactions in the disk 
give rise to an effective viscosity and allows the transport 
of angular momentum and mass in the gas disk. In our 
scenario turbulence is assumed to be generated by insta- 
bilities involving self-gravitation and to be maintained by 
the energy input, which is provided by differential rota- 
tion of the disk and mass transfer to smaller galactic radii 
via cloud-cloud interactions. 

A description of the turbulent viscosity is given, which 
is formally equivalent to a prescription derived from the 
collisional Boltzman equation. The vertical gravitational 
acceleration in the gas disk is either due to a (i) domi- 
nating central mass, (ii) dominating stellar disk mass, or 
(iii) the vertically self-gravitating gas disk itself. For all 
three cases we derive analytical expressions for disk pa- 
rameters (density, surface density, velocity dispersion, disk 
height, driving wavelength, and viscosity) as functions of 
the Toomre parameter Q, the turbulent Reynolds number 
Re, the mass accretion rate M, the galactic radius R, and 
the angular velocity O of the gas. 



The model does not include radiation. The "thermo- 
stat" equation of the standard accretion disk model is re- 
placed by an equation which assumes that the energy flux 
generated by mass inflow and differential rotation is en- 
tirely transported by turbulence to smaller scales where 
it is dissipated. Whereas case (iii) corresponds to the a 
Q ^ 1 disk, case (i) and (ii) require Q > 1. 

The structure of the resulting clumpy gas disks allows 
us to derive global volume filling factors of self-gravitating 
clouds as a function of Q and Re. Along the same line 
analytical expressions for the molecular fraction and the 
star formation rate are given. 

On the basis of our analytical model we conclude that 

1. constant velocity dispersions in the gas as observed for 
spiral galaxies can be reproduced by models with a (i) 
dominating central mass and (iii) a vertically selfgrav- 
itating gas disk. 

2. In models (ii) with a dominating stellar disk the veloc- 
ity dispersion depends on the rotation velocity. 

3. The driving wavelength of the turbulence equals ap- 
proximately the disk scale height if the vertical accel- 
eration is not dominated by the stellar disk. 

4. The effective turbulent viscosity v depends on the 
galactic radius. For a disk in Keplerian rotation with 
a (i) dominating central mass v oc R^. The other two 
model classes (ii) and (iii) give i/ oc i? for a constant 
velocity rotation curves. 

5. The molecular gas surface density and the star forma- 
tion rate are proportional to the square of the angular 
velocity Sjja , oc in self-gravitating gas disks. 

The model of a selfgravitating disk in vertical direction 
is applied to the Galaxy and gives a good description of 
its gas disk. Furthermore, we investigate the role of the 
galaxy mass for the morphological classification of spiral 
galaxies. 

A comparison of our analytical model to spiral galaxies 
shows that 

1. the observed global star formation time scale implies 
a Reynolds number Re ^ 50. 

2. our model is consistent with observations. That is 
the model reproduces the physical parameters of the 
Galaxy as derived from observations. 

3. the derived mass accretion rate of the galaxy is M ~ 
10~^ MQyr~^. It is thus much smaller than the star 
formation rate M*/M ~ 100. 

4. the dependence of several physical parameters of spiral 
galaxies on morphological type might be at least partly 
due to the mass-morphology relation, i.e. galaxies of 
later types are more massive. 

In our model we assume that turbulence is generated by 
local instabilities due to selfgravitation. The energy input 
from accretion and differential rotation is high enough to 
maintain the ISM turbulence. Moreover, the energy dis- 
sipation rate per unit mass due to the self-gravitation of 
the gaseous disk in z direction has the same order of mag- 



14 



B. VoUmcr and T. Bcckert: Turbulent viscosity in clumpy accretion disks 



nitude as the energy per unit time and unit mass trans- 
ported to smaller scales by turbulence. Thus, in terms of 
energy conservation, turbulence generated and maintained 
by gravitation can account for the viscosity of spiral galax- 
ies. This does not exclude energy input in the turbulence 
through supcriiovae. Nevertheless, our model where accre- 
tion in the galactic potential is sufficient to maintain ISM 
turbulence is fully consistent with observations. 
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